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Abstract 

We  investigated  the  transmission  of  Dirac  electrons  through  a  potential  barrier  in  the  presence 
of  circularly  polarized  light.  An  anomalous  photon-assisted  enhanced  transmission  is  predicted 
and  explained.  It  is  demonstrated  that  the  perfect  transmission  for  nearly  head-on  collision 
in  infinite  graphene  is  suppressed  in  gapped  dressed  states  of  electrons,  which  is  further 
accompanied  by  a  shift  of  peaks  as  a  function  of  the  incident  angle  away  from  head-on  collision. 

In  addition,  the  perfect  transmission  is  partially  suppressed  by  a  photon-induced  gap  in  illumi¬ 
nated  graphene.  After  the  effect  of  rough  edges  of  the  potential  barrier  or  impurity  scattering 
is  included,  the  perfect  transmission  with  no  potential  barrier  becomes  completely  suppressed 
and  the  energy  range  for  the  photon-assisted  transmission  is  reduced  at  the  same  time. 

(Some  figures  may  appear  in  colour  only  in  the  online  journal) 


1.  Introduction 

Research  activity  on  the  properties  of  graphene  has  been 
growing  rapidly  ever  since  its  experimental  discovery  and 
the  demonstration  of  its  unusual  properties  arising  from  its 
energy  band  structure  [1,2].  The  novel  properties  of  graphene 
may  be  attributed  to  its  massless  Dirac  fermions  at  the  Fermi 
energy  [3],  An  interesting  consequence  of  the  Dirac  electron 
is  the  Klein  paradox  [4],  in  which  an  electron  in  graphene 
undergoes  unimpeded  tunneling  through  potential  barriers  of 
arbitrary  height  and  thickness.  This  property  of  the  Dirac 
electrons  is  due  to  their  linear  energy  dispersion  relation  or 
helicity.  Electrons  are  said  to  be  chiral  if  their  wavefunctions 
are  eigenstates  of  the  chirality  operator  h  =  a  •  p/(2 p),  where 
o  =  {<yx,  erv}  is  the  Pauli  vector  consisting  of  Pauli  matrices 
and  p  =  {px,  py]  is  the  electron  momentum  in  he  graphene 
layers.  Electrons  in  graphene  near  the  K  points  (around  the 
corners  of  the  hexagonal  Brillouin  zone)  are  chiral  due  to  the 
fact  that  the  chirality  operator  is  proportional  to  the  Dirac 
Hamiltonian  which  automatically  makes  chirality  a  good 
quantum  number. 

The  fact  that  the  Klein  paradox  is  also  obtained  in  bilayer 
graphene  makes  this  effect  even  more  sophisticated.  This 
leads  us  to  realize  that  the  Klein  paradox  is  not  simply  due  to 
linear  electron  dispersion  but  may  be  observed  for  both  mass¬ 


less  and  massive  quasiparticles  [4],  In  this  paper,  we  consider 
a  sharp  p-n-p  junction  or  potential  barrier  profile.  This  type  of 
potential  can  be  constructed  by  an  underlying  metal  contact  or 
an  insulating  strip  [5]  and  has  been  employed  to  demonstrate 
unimpeded  transmission  [3],  Apart  from  using  an  optical  field, 
one  may  also  create  potential  barriers  (p-n  junctions)  as  well 
as  energy  gaps  and  coherence  using  the  surface  curvature  as 
discussed  in  the  paper  by  Atanasov  and  Saxena  [6], 

Both  the  tight-binding  model  and  the  k-p  approximations 
for  an  infinite  graphene  sheet  accurately  show  that  electrons 
and  holes  have  linear  energy  dispersions  s(k)  =  shvpk  = 

shvpyjkx  +  ky  with  no  gap,  where  s  is  the  electron-hole  parity 
with  s  =  1  for  electrons  and  s  =  —  1  for  holes.  For  the  potential 
barrier  described  above,  there  is  a  translational  symmetry  in 
the  y  direction  parallel  to  the  boundaries  of  the  potential  so 
that  ky  is  conserved.  In  contrast,  the  longitudinal  component 
kx  is  modified  by  the  potential  so  that  when  the  particle 

has  energy  E  we  have  kxj  =  (/zi’p)-1  J~(E  —  Vq)2  —  (hv^ky)2, 
where  V,  is  the  potential  in  the  region  i. 

The  electron  effective  mass  and  its  properties  will  play 
a  crucial  role  in  our  analysis.  In  infinite  intrinsic  graphene, 
the  Dirac  electron  is  massless  although  the  mass  can  be 
introduced  [7]  or  implemented  experimentally.  For  bilayer 
graphene,  the  particle  effective  mass  exists  in  any  possible 
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approximation  [8],  leading  to  the  existence  of  evanescent 
terms  in  the  wavefunction  [4],  although  this  does  not  violate 
chirality  symmetry  and  the  Klein  paradox. 

The  outline  of  the  rest  of  this  paper  is  as  follows. 
In  section  2,  we  derive  the  eigenstates  and  transmission 
coefficient  for  Dirac  electrons  dressed  with  photons  from  a 
circularly  polarized  light.  Section  3  is  devoted  to  a  formalism 
in  which  the  roughness  of  the  boundaries  for  the  potential 
barrier  is  included  phenomenologically  in  calculating  the 
transmission  coefficient;  numerical  results  are  also  presented. 
A  brief  summary  is  given  in  section  4. 

2.  Transmission  coefficient  for  dressed  electron 
states 


and  the  Jaynes-Cummings  Hamiltonian 

„  +  „  w  „ 

'H\-r  =  h(ona  a - —  (cr+a  +  er_a  ),  (5) 

2  VN 

which  corresponds  to  a  two-level  quantum  optical  system 
and,  most  importantly,  can  be  solved  analytically.  Here,  we 
have  defined  k±  =  kx  ±  i ky  and  N  represents  the  number 
of  radiation  quanta  (intensity)  for  the  incident  optical  field. 
We  only  consider  the  situation  such  that  the  electron-photon 
interaction  amplitude  vv  is  much  less  than  either  the  photon  or 
the  Dirac  electron  energy,  i.e., 


/  2Nhe2vl 

w  —  2.  - —  =  2a  two  <  h(0Q,  (6) 

V  e0tt>oV 


It  was  shown  recently  [9,  10]  that  when  Dirac  electrons  in  a 
single  graphene  layer  are  interacting  with  an  intense  circularly 
polarized  light,  the  electron  states  will  be  dressed  by  photons. 
The  main  idea  of  the  present  study  is  to  investigate  the 
transmission  properties  of  such  dressed  electrons  for  the  case 
of  single  layer  graphene.  We  go  beyond  the  approximations 
used  in  [9]  by  retaining  the  results  up  to  the  order  of 
(D( A4)  so  that  we  are  able  to  investigate  the  difference 
between  the  dressed  states  and  the  massive  Dirac  electrons 
described  below  by  the  Hamiltonian  in  (24).  Here,  A  is  a 
quantity  measuring  the  induced  gap  between  the  valence  and 
conduction  bands  of  the  dressed  electrons. 

We  begin  with  the  electron-photon  interaction  Hamilto¬ 
nian 

U  =  vfct  •  (p  -  eAcirc),  (1) 

where  vp  is  the  Fermi  velocity  and  the  vector  potential  for 
circularly  polarized  light  of  frequency  coq  can  be  expressed  as 


where  a  =  w/hcoo  1.  The  two  eigenstates  of  the 
Hamiltonian  (5)  could  be  obtained  with  an  expansion  over  the 
basis  of  just  two  functions  1 1}^  =  |  -f ,  N)  and  1 2)N  =  |  j, ,  N+ 1 ) 
for  each  N  value,  that  is, 

I'l'tw)  =  Hn\1)n  +  Dv|2);v,  (7) 

=  Hn\2)n  -  VArl  1)jv,  (8) 

which  corresponds  to  the  following  Hamiltonian  in  the  chosen 
basis  |  [)N  and  |2)jy: 


/qi  fii2 

N 

_a  /AH-  1 

=  hw  0 

T2  V  N 

1*21  1*22 

a 

LT2 

/AH-  1 

N+  1 

V  N 

In  this  way,  the  transformation  (7)  and  (8)  becomes  just  a 
simple  rotation  in  the  Hilbert  space.  This  expansion  yields  the 
eigenvalue  equation  and  its  solution  below: 


Arirc  - 


eo  coqV 


(e+a  +  e_fi^) 


h 


2eo  a>0y 


[(«  +  flT)eA-  +  i (a  -  aT)ey] 


(2) 


in  terms  of  photon  creation  and  destruction  operators  a1 
and  afi ,  respectively.  Here,  V  is  the  mode  volume  of  an 
optical  field.  In  order  to  study  the  complete  electron-photon 
interacting  system,  we  must  add  the  field  energy  term  hcooa^a 
to  the  Hamiltonian  (1).  As  usual,  we  seek  the  wavefunction  in 
the  form  of  a  plane  wave  'T(r)  =  elk'ri//(fi),  which  results  in 
the  following  reduced  Hamiltonian: 


=  tuDo  o'  a  +  fivpcr  •  k  — 


2/ie2Vp 

Cr\ 


(cr+a  +  cr-fi1),  (3) 


where  cr+  =  \(ax  +  ierv)  =  and  a_  =  \(ax  -  i ay)  = 

^  .  The  reduced  Hamiltonian  (3)  for  infinite  graphene  can 

be  understood  as  the  sum  of  two  parts,  namely,  the  Dirac 
Hamiltonian 


Dirac  =  frv F  CT  •  k  =  tlVF(crxkx  +  Gyky) 

=  /ivp((T_k_|_  +  cr+k_) 


(a  lN+  1 V 

(e  -  Nhw 0)[e  -  (AH-  1  )tico0\  -  I  -  hco0  y  ^  )  =  0, 

(10) 

=  (N + 0 r,m  *  T!\/l+“2(^) 

(  1  1  \  a2 

-  I  A^+  -=F  ^  —Hcoq,  (ID 

where  e-|-  and  ej,  correspond  to  the  lower  +  and  upper  — 
signs  in  the  solution.  For  simplicity,  we  assume  here  that  the 
radiation  is  classically  strong  with  N  1.  However,  we  note 
that  we  may  set  N  +  1  ±2  N  only  in  the  terms  0(a2)  but  not 
in  the  terms  with  Nhm 0.  Furthermore,  it  is  a  simple  matter  to 
obtain  the  expansion  coefficients  /z,#  and  %  as 


Hn  —  cos  9C  —  1 - , 


tan0r  = 


vN  =  sin  0C  - 

aJ(N+  \)/N 
1  +  y/\  +a2(N+  1  )/N 


(12) 


One  can  easily  verify  that  the  wavefunctions  (7)  and  (8)  are 
the  eigenstates  of  the  Jaynes-Cummings  Hamiltonian  (5)  with 


(4) 
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energies  (11),  to  any  order  of  a.  Here,  for  all  the  above 
derivations,  we  have  employed  the  standard  relations,  i.e., 

«t|t  i,N)  =  s/N+  l|t  \,N+  1), 
h\\  \,N)=Vn\\  \,N-\), 
a+U,N)  =  \f,N),  <t+\\,N)  =  0, 

=  U,N),  a-\i,N)  =  0. 

We  now  look  for  the  eigenstates  of  the  Hamiltonian  (3)  as 
expansions  over  the  set  of  Jaynes-Cummings  Hamiltonian 
eigenfunctions  (7)  and  (8).  We  will  confine  our  attention  to 
the  field  source  with  only  three  nearest  photon  occupation 
numbers,  i.e.,  N  —  Nq  —  l.  No  and  No  +  1,  leading  to 


Wo+l 

|4>W>-  £  (Cu(k)  |'Ptl/)+C2,/(fc)|^>).  (14) 
e=N0-i 

We  know  that  the  eigenfunctions  (7),  (8)  corresponding  to 
different  numbers  N  are  orthogonal  to  each  other.  First  acting 
the  Hamiltonian  (3)  on  (14),  and  then  multiplying  both  sides 
of  the  expansion  (14)  by  ('I'-f,At0|,  ('I'tWo+il  and 

(4'|,v0l,  this  results  in  the  following  four  equations: 


Ci, N0  (M)  - 
=  eCi,iv0, 


or  Nq  +  1 


C2,n0-i  No  + 


No 

r  No 


hcoo  +  /tvF('I'|jVolcr  •  k|  OfA')) 


hcoo 


C\,N0+i 


4  N0  -  1 , 

+  fiVF('k|iA?0_i  I CT  •  k|  <$>(A'))  =  eC2.At0-l, 
a2  No  +  2" 


(No  +  1)  - 


hco  o 


C2, No 


4  Nq  +  1  _ 

+  /fVF^tJVo+il*7  •  k|<t>(k))  =  eCi,At0+i, 
a2  No  +  1  " 


(No  +  1)  + 


hco  o 


4  No 

/jvF('k^iv0|a  •  k|d>(k)>  =  sC2,n0- 


(15) 


For  chosen  number  No,  each  pair  of  these  equations  describes 
two  energy  subbands  that  are  separated  by  ±2 a2  hcoo /2  at  k  = 
0.  Here,  we  include  only  the  four  nearest  energy  subbands, 
corresponding  to  the  dressed  states  with  different  photon 
occupation  numbers  N  and  electron  states  with  subband 
indices  •f  (lower  energy)  and  J,  (higher  energy). 

Taking  into  account  the  following  relations  for  the  Dirac 
Hamiltonian: 


ct  •k|T'tJv0>  =  fJ.Nok+\i,N0)  -  vjv0k_ 1 1 , IVo  +  1).  (16) 
ct  -k|^0)  =  -/x;v0k_|t,77o  +  1)  -  UVo^+U-M)}, 

and  with  the  simplifications  described  above,  we  can 
explicitly  write  out  the  Dirac  Hamiltonian  terms  in  equations 
(15)  and  (15).  The  presence  of  photon  occupation  numbers 
does  not  follow  from  the  dressed  states  Hamiltonian  and 
cannot  be  determined  just  from  the  ratio  between  the 
photon  field  energy  Nhcoo  and  the  electron-photon  interaction 
amplitude  w.  It  should  result  from  the  model  of  the  circularly 
polarized  light  source.  If  we  consider  electrons  near  the  K 
point  such  that  the  kinetic  energy  hvp  k  ri  a2hco <j/2,  we  can 
retain  only  three  photon  occupation  numbers  Nq  —  l,  Nq  and 


Nq  +  1 ,  which  is  consistent  with  the  approximation  described 
above.  Consequently,  we  arrive  at  a  system  which  consists 
of  weakly  coupled  equations  to  determine  the  two  nearest 
subbands,  i.e., 


Noha>o  —  A 

—  p2frl'pk_ 

/rvfivpk+ 

0 

fj,2hvpk+ 

NqHcoq  +  A 

0 

—v/ohvpk+ 

—v/jihv^k- 

0 

(No  +  l)fift>o  - 

A  —fi^hvpk- 

0 

v/r/7Vpk_ 

/r2/n’pk+ 

(No  +  1  )ha>o  +  A 

"Cf 

"Cf 

c2 

c2 

C3 

—  £ 

C3 

_C4_ 

c4_ 

(17) 


where  fi  =  /rjv0,  v  =  VfjQ,  A  a2hcoo/4,  C\  =  Ci,jv0,  C2  = 
C2,No—ii  C3  —  Ci,jv0+i  and  C4  =  C2,n0-  This  leads  to  the 
following  energy  dispersions  with  /r  ss  1 : 


£l,l(k)  =  (Nq  +  j)k(Do 

-  \  hco o^l  +  r\  +  | (1  +  v2)k2  ±  2^(f  v2k2  +  1  )(^k2  +  17), 


(18) 


£3,4(&)  =  (Nq  +  \)h(D  0 


+ 


;  heoo^l  +  1?  +  |(1  +  v2)k2  T  2  J(^v2k2  +  1  )(fk2  +  17), 


(19) 


where  £  =  ( 2vp/a>o )2  and  11  —  (2A/ha>o)2  ~  a4/4. 

Another  approximation  which  we  may  used  here  is  to 
consider  only  small  wave  vectors  k  so  that  we  get  two 
independent  pairs  of  subbands  separated  by  energy  hco 0.  Let 
us  estimate  the  order  of  magnitude  for  each  off-diagonal 
element  in  the  Hamiltonian  (17).  Using  the  relations  in 
equations  (12),  we  know  that  the  expansion  coefficients  are 
li  ^  1  and  v  ^  a,  where  «  1  is  the  ratio  between  the 

electron-photon  interaction  energy  and  the  photon  energy 
hcoo ■  Note  that  the  energy  gap  is  A  =  a2  hcoo/ 4  ^  0(a2).  Now 
let  us  consider  small  momenta  k  such  that  hvpk  orhcoo. 
As  a  result,  we  find  that  ±fi2hvpk±  ^  orhcoo  ^  0(a2)  and 
±(ivhvpk±  a2 hcoo  ^  0(a2)  — >  0.  After  all  the  0(a 3)  terms 
in  equation  (17)  are  neglected,  we  obtain 


Nohcoo  —  A 

— /x2Avpfe_ 

0 

0 

/r2/ivpk+ 

Nohcoo  +  A 

0 

0 

0 

0 

(No  +  l)ha>o  — 

A  — /x2fjvpk_ 

0 

0 

fi2hvpk+ 

(No  +  l)ha>o  +  A 

"Cf 

"cf 

C2 

c2 

C3 

—  £ 

C3 

_C4_ 

c4_ 

(20) 


This  gives  us  two  independent  pairs  of  energy  subbands 
and  the  corresponding  two  independent  parts  of  the 
wavefunction.  Each  independent  part  is  a  two-component 
spinor  wavefunction,  where  each  component  consists  of  two 
independent  terms. 
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The  corresponding  eigenvalue  equation  for  the  case  of 
two  independent  pairs  of  subbands  yields,  as  expected, 

{hr v£  Ic  +  A2  -  (e  -  No  hco o)2}  {h24  Ic  +  A2 

-  [e  -  (Nq  +  1)  hcoo]2}  =  0.  (21) 

If  we  consider  the  transmission  of  such  states  through  a 
potential  barrier  whose  height  is  equal  to  hco o,  one  of  the 
terms  in  region  ‘1’  will  exactly  match  the  other  in  region 
‘2’  (the  potential  region)  so  that  a  part  of  the  wavefunction 
will  be  completely  transmitted.  This  will  lead  to  a  substantial 
increase  in  the  total  transmission.  On  the  other  hand,  if  we 
consider  a  larger  number  of  such  subband  pairs,  this  effect 
will  not  be  significant  since  this  complete  transmission  occurs 
only  for  one  of  the  wavefunction  terms.  This  is  reminiscent 
of  the  Klein  paradox  since  it  does  not  depend  on  the  barrier 
width  as  long  as  the  barrier  height  is  exactly  equal  to  hco o- 
The  revelation  of  such  an  anomalous  increase  for  the  dressed 
state  tunneling  is  an  important  discovery  of  the  present  paper. 

The  simplest  possible  approximation  is  to  keep  only  two 
terms  with  the  coefficients  [C\(k),  C2(k)]  leading  to  the  two 
nearest  energy  subbands  separated  by  a  gap  2  A  ±2  a2  hcoo /2  at 
k  =  0  due  to  electron-photon  interaction.  The  block-diagonal 
approximation  introduced  in  equation  (20)  is  valid  only  for  a 
single  value  Nq  of  photon  number  as  well  as  small  values  of 
k.  In  this  situation,  the  other  subband  pairs  become  far  away 
from  the  considered  one  and,  therefore,  their  influence  can  be 
neglected.  We  will  use  this  approximation  in  the  rest  of  our 
paper. 

Equation  (20)  results  in  a  simplified  algebraic  system 
determining  the  pair  of  coefficients  C\ ( k )  and  C2(k),  i.e., 

(Nohcoo  ~  A)  Ci  +  hv-p(kx  +  iky)C2  =  £  C\ , 

(. Nohcoo  +  A)  C2  +  tivp  (kx  —  iky)  C\  =  eC2. 

The  non-trivial  solution  of  these  equations  gives  the  energy 
dispersion  which  was  previously  obtained  in  [9]  as 

s(k)  —  Nq  hco o  ±  J  A2  +  h2v^.k2  =  No  hcoo 

+  2  +  h2v2k2.  (23) 

If  the  electron-photon  interaction  is  removed,  then  a  = 
w/2  ->  0,  and  the  energy  dispersion  relations  (23)  demon¬ 
strate  a  non-interacting  system  consisting  of  a  Dirac 
electron  £p(k)  =  fihv f  \k\  and  photons  Nohcoo .  All  the  other 
energy  subbands  are  separated  at  least  by  hcoo  'T>  w  and 
can  be  neglected,  which  justifies  the  above  two-subband 
approximation.  In  this  notation,  =  ±1  is  the  dressed 
conduction/valence  band  index  corresponding  to  bare  elec¬ 
tron/hole  bands  for  infinite  graphene  when  A  — >  0. 

The  system  (22)  is  formally  similar  to  the  eigenvalue 
equations  for  the  case  of  the  effective  mass  Dirac  Hamiltonian 


TL  —  hyp  a  •  k  +  V(x) 


1  0 
0  1 


+  A  cr3, 


(24) 


where  03  is  a  Pauli  matrix  and  V(x)  is  a  one-dimensional 
potential.  The  electron  dispersion  and  transmission  properties 
for  both  a  single  as  well  as  multiple  square  potential 


barriers  have  been  studied  [7,  11]  for  monolayer  and  bilayer 
graphene  [12].  It  was  also  shown  that  a  one-dimensional 
periodic  array  of  potential  barriers  leads  to  multiple  Dirac 
points  [13].  Several  papers  have  introduced  an  effective  mass 
term  into  the  Dirac  Hamiltonian  for  infinite  graphene  which 
may  be  justified  based  on  different  physical  reasons  [14]. 
For  example,  it  has  been  shown  [15]  that  an  energy  bandgap 
in  graphene  can  be  created  by  a  boron  nitride  substrate, 
resulting  in  a  finite  electron  effective  mass.  However,  we 
emphasize  that  the  analogy  between  the  Hamiltonian  (24)  and 
that  for  irradiated  graphene  is  not  complete  since  that  would 
correspond  to  A  <  0.  Although  this  difference  does  not  result 
in  any  modification  of  the  energy  dispersion  term  containing 
A2,  it  certainly  modifies  the  corresponding  wavefunction.  The 
occurrence  of  a  tunable  energy  gap  in  graphene  coupled  to  a 
laser  source  is  discussed  in  [16]. 

It  was  demonstrated  in  [17,  18]  that  strained  superlattices 
can  be  used  to  open  an  appreciable  energy  gap  in  the  graphene 
energy  spectrum  and  modify  the  spin-orbit  interaction 
considerably.  However  the  advantage  of  using  a  circularly 
polarized  light  is  that  the  energy  gap  is  tunable  by  varying 
the  frequency  and  intensity  of  the  applied  optical  field.  In 
addition,  we  know  that  an  energy  gap  can  also  be  created 
by  using  the  surface  curvature  [6].  However,  for  the  case 
of  curved  or  strained  graphene,  the  longitudinal  momentum 
is  no  longer  a  good  quantum  number  and  the  translational 
symmetry  is  lost.  As  a  result,  the  formalism  presented  in  this 
paper  for  finding  the  transmission  amplitude  can  no  longer  be 
applied  to  this  case. 

Instead,  we  consider  another  case  in  which  the  energy 
gap  is  created  by  a  substrate  and  a  circularly  polarized  light 
at  the  same  time.  In  this  case,  we  obtain  the  following 
Hamiltonian  [19]: 


hcoo  »  „+ 

TL  =  hcoo  a'  a  +  hv f  a  •  k  -| — —  <73  —  ahcoo(a+a  +  o-a  ) 


„  +  A  hcoo 

=  hcoo  a  a  +  HDirac  H — —  er3  +  7/j-c, 


(25) 


which  gives  rise  to  the  energy  dispersion 


e(k)  =  y/(hvFk)2+  A2ub  + A2.  (26) 

Here,  Asub  is  the  substrate-induced  energy  gap  and  A  is  the 
energy  gap  induced  by  the  electron-photon  interaction.  The 

total  energy  gap  is  then  given  by  Atotai  =  ^/^sub  +  where 
we  only  consider  the  case  with  two  nearest  energy  subbands. 

The  interaction  between  the  Dirac  electrons  in  graphene 
and  a  circularly  polarized  light  has  been  considered  in  the 
classical  limit  in  [20].  In  this  limit,  a  gap  in  the  Dirac 
cone  opens  up  due  to  nonlinear  effects.  The  dressed-state 
wavefunction  has  the  form 


4>dr(£)  = 


Ci(k) 
pC2(k)e * 


with  C\(k)  ^  C2(k)  given  by 


Cf(k) 


1 

V/2(l  +  K2)T2yVl  +  y2 


(27) 
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C9  ( k )  =  zb 


yr+ 


y  t  y 


j2(l  +  Y2)T2yVi  +  y2 


(29) 


corresponding  to  the  energy  subbands  s(k )  =  Noticoo  ± 
J  A2  +  k2.  In  this  notation,  y  —  A  /(/H’f  A:)  and  0  is  the 

angle  which  k  makes  with  the  longitudinal  x  axis.  Without  an 
optical  field,  i.e.  A  =  0,  we  obtain  Cf  —  =  1/72.  In  the 

limit  (A  <5jC  k),  the  coefficients  exhibit  peculiar  symmetry  with 


y  y2  3y 3 

2V2  8V2  T  I6V2’ 

y  r_  +  3  y3 

2V2  T  8V2  1672' 


(30) 


We  note  that  this  expansion  is  not  valid  too  close  to  the  Dirac 
point  and  should  not  be  used  for  an  arbitrary  wave  vector  to 
calculate,  for  example,  the  polarization  function.  Additionally, 
one  may  verify  that  C\(k)  ^  C2(k)  for  any  chosen  A  in  the 
range  of  validity.  Consequently,  the  chiral  symmetry  is  broken 
for  electron  dressed  states 


h^drik)  = 


1  a  ■  p 

C\(k) 

1 

PC2(k) 

2  P 

PC2(k)  e'* 

“  2 

Ci  (it)  e1^ 

Clearly,  it  follows  from  (31)  that  the  non-chirality  of 
the  dressed  electron  states  becomes  significant  if  the 
electron-photon  interaction  (the  leading  y  term)  is  increased. 
This  affects  the  electron  tunneling  and  transport  properties. 
We  now  turn  to  an  investigation  of  the  transmission  of  electron 
states  through  a  potential  barrier  when  graphene  is  irradiated 
with  a  circularly  polarized  light. 

For  simplicity,  we  consider  a  square  potential  barrier  of 
height  Vo  given  by  V(x)  =  Vo[9(x)  —  9(x  —  Wo)],  where  Wo 
is  the  barrier  width  and  6(x )  is  the  Heaviside  step  function. 
Since  the  wave  number  kx  is  the  same  in  regions  1  (x  <  0)  and 
3  (x  >  Wo)  and  the  current  component  is  jx  =  O :  <rA<t>,  we  only 
need  the  wavefunction  continuity  at  the  potential  boundaries 
for  the  system  considered.  From  this  continuity  condition,  the 
transmission  probability  T  can  be  determined  from  T  =  |f|2, 
where  t  is  the  transmission  coefficient  or  the  amplitude  of  the 
wave  propagating  forward  in  region  3.  Here,  we  only  show 
an  analytical  expression  for  the  transmission  coefficient  of 
the  dressed  states  in  the  limit  of  £  «  ko  corresponding  to 
kx\  kx 2: 

cos  20 

T  =  - — - 

cos2(kA-2Wo)cos20  +  sin2  (£*2  Wo) 

_ 3cos2 (kx 2  Wq)  cos20sin2 (kx2  Wp) 

Az2 Vpk2i [sin2 (kx2 Wo)  +  cos 2(kx2W0)  cos20] 

Here,  6  =  tan-1  (ky/kx2)  — >•  0,  0  =  tan-1  (ky/kx\),  the  second 
term  includes  the  effect  of  electron-photon  interaction  (ocA2) 
to  the  leading  order.  In  addition,  we  only  show  those  terms  of 
the  lowest  order  in  e/Vo-  There  is  another  relevant  study  [21], 
investigating  the  tunneling  of  Dirac  electrons  with  a  finite 
effective  mass,  a  parabolic  dispersion  in  the  presence  of 
an  energy  gap,  and  a  certain  chirality,  through  a  potential 
region.  In  that  study,  a  particle  tunneling  through  a  square 


potential  barrier  differs  from  both  the  Dirac  electrons  and  the 
dressed  states  of  electrons  under  a  circularly  polarized  light 
illumination. 

For  nearly  head-on  collision  with  ky  kx  1  kx2  for 
high  potential  as  well  as  for  infinite  graphene  (A  — »■  0),  the 
transmission  coefficient  has  the  following  simplified  form: 

T  =  1  —  sin2 (kx2Wo)(d2  -  2 pd<p  +  02),  (33) 

where  we  assume  Vo  ^  s,  0  0  <£;  1  and  /3  =  ±1. 

3.  Numerical  results  and  discussion 

In  our  numerical  calculations,  energies  will  be  measured 
in  units  of  (3kfCit/2)  with  the  carbon-carbon  distance  a  & 
1.42  A  and  the  hopping  parameter  t  =  2.773/2  eV.  We 
measure  the  wave  vector  in  units  of  the  Fermi  wave  number 
kp  and  write  its  components  as  kx  =  cos  0  and  ky  =  sin  0  in 
terms  of  the  angle  of  incidence  0. 

In  figure  1,  we  present  the  transmission  for  dressed 
electron  states  with  arbitrary  energy  and  angle  of  incidence. 
We  clearly  see  that  dressing  ruins  the  Klein  paradox  in 
(a)  for  head-on  collision  with  0  =  0.  The  resonant  peaks 
are  shifted  for  the  other  incoming  angles  in  (b)  and  the 
effect  is  stronger  for  small  incident  angles.  In  (c)  and  (d), 
the  transmission  probability  plots  are  given  in  terms  of  the 
longitudinal  momentum  kx\  in  front  of  the  barrier  and  kx2  in 
the  barrier  region.  We  find  that  the  intensities  and  locations 
of  the  transmission  peaks  in  (d)  are  distorted  compared  to 
those  for  infinite  graphene  in  (c).  The  diagonal  kx\  =  kx2 
corresponds  to  the  absence  of  a  potential  barrier  and  should 
yield  a  complete  transmission  for  A  =  0.  However,  the 
condition  (-/ (e  —  Vo)2  —  A2  >  hv f  ky)  for  dressed  states  must 
be  satisfied,  which  makes  the  diagonal  incomplete  (missing 
diagonal  for  small  kx  \  and  kx2 )  due  to  the  occurrence  of  an 
induced  gap. 

Figure  2  displays  the  effect  due  to  the  electron-photon 
interaction  on  the  electron  transmission  in  terms  of  incoming 
particle  energy  s  and  angle  of  incidence  0.  From  the  figure, 
we  see  the  Klein  paradox  as  well  as  resonant  tunneling  peaks 
in  the  transmission  probability  for  regular  infinite  graphene 
with  A  =  0  in  (a).  The  dark  ‘pockets’  on  both  sides  of  f:  —  Vo 
in  (a)  demonstrate  zero  transmission  for  the  case  |e  —  Vo|  <$( 
e,  which  results  in  imaginary  longitudinal  momenta  kx2  for 
most  incident  angles  and  produces  a  completely  attenuating 
wavefunction.  When  a  small  gap  is  opened  in  (b)  for  the 
dressed  states  of  electrons  in  graphene  under  illumination  by 
a  circularly  polarized  light,  we  observe  a  set  of  complete 
transmission  branches,  where  a  strong  dependence  on  0  for 
the  lower  branches  is  seen.  However,  this  0  dependence  is 
greatly  suppressed  when  the  dressed-state  gap  is  increased 
in  (c),  leaving  us  a  set  of  equally  distant  branches  due  to 
photo-assisted  electron  tunneling. 

As  mentioned  above,  by  including  more  than  the  two 
nearest  subbands,  as  shown  in  figure  3(a),  the  electron 
dispersion,  wavefunction  and  transmission  amplitude  will 
be  modified.  According  to  the  approximation  adopted  by 
Kibis  [9],  the  two  subband  pairs  may  be  considered 
as  independent.  Under  the  condition  of  equal  transverse 
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Figure  1.  Transmission  probability  T  of  electrons  as  functions  of  incoming  energy  e,  incident  angle  <p,  and  longitudinal  wave  numbers 
hv p  kx i  and  /ivp  kx.  In  (a),  the  transmission  probability  (purple  dashed  curve)  is  plotted  as  a  function  of  the  incident  particle  energy  for 
head-on  collision  (<j)  =  0).  The  red  solid  line  is  the  transmission  probability  T  =  1  (Klein  paradox)  for  infinite  graphene  when  the  energy 
gap  A  =  0.  In  (b),  we  show  the  angular  distributions  of  T  for  infinite  graphene  (solid  red  curve)  and  the  dressed  states  (purple  dashed  curve) 
with  e  =  Fo/6.  In  (c)  and  (d),  the  transmission  probability  is  plotted  as  a  function  of  the  longitudinal  momenta  in  regions  ‘1’  and  ‘2’  (before 
the  potential  barrier  and  in  the  barrier  region,  respectively)  for  infinite  graphene  (c)  and  irradiated  graphene  (d).  The  diagonal  kx\  =  kx 2 
corresponds  to  the  absence  of  a  potential  barrier  (Vo  —  0)  and  displays  a  perfect  transmission  with  T  =  1  as  long  as  both  longitudinal 
momenta  are  real  determined  by  both  the  barrier  height  Vo  and  the  energy  gap  A. 


Figure  2.  Transmission  probability  T  as  functions  of  both  the  incoming  electron  energy  e  and  the  angle  of  incidence  <j>.  Plot  (a)  is  for 
infinite  graphene  (A  =  0)  with  the  obvious  Klein  paradox  for  0  =  0.  Plots  (b)  and  (c)  show  the  transmission  T  for  electron  dressed  states,  in 
the  two  nearest  subbands  approximation,  with  gap  energies  A  =  Vo/15  (b)  and  A  =  Vo/5  (c). 


momentum  ky  for  both  terms  of  the  particle  wavefunction, 
the  first  term  of  the  wavefunction  in  region  ‘2’  with  a 
potential  barrier  is  similar  to  the  second  term  in  the  regions 
without  potential.  Therefore,  the  states  exactly  match  across 
the  potential  boundary,  which  should  definitely  increase  the 
total  transmission  amplitude.  The  other  possibility  is  that 
incoming  angle  and  momentum  of  the  second  term  are  totally 
independent  of  the  first  one  and  lead  to  resonant  transmission 
regardless  of  the  transmission  amplitude  of  the  first  term. 
Based  on  our  derived  results,  we  find  that  the  transmission 
should  increase  even  for  non-split  energy  subband  pairs 
and  corresponding  wavefunctions.  Since  the  existence  of 
certain  photon  occupation  numbers  N  is  determined  by  the 
laser  source,  we  can  consider  that  only  one  pair  of  the 


dressed  state  subbands  is  occupied  while  the  other  subband 
pairs  are  unpopulated.  As  illustrated  in  figure  3(b),  in  the 
potential  region  ‘ 2’,  a  particle  may  populate  another  subband 
corresponding  to  a  different  number  of  photons  instead 
of  changing  its  longitudinal  momentum  for  barrier  heights 
exactly  equal  to  a  multiple  of  Hcoq.  The  opposite  transition 
will  occur  at  the  boundary  between  regions  ‘2’  and  ‘3’.  This 
will  result  in  unimpeded  tunneling  T  =  1  independent  of  the 
barrier  width,  as  seen  from  figure  3(c),  which  is  expected  to 
be  a  major  contribution  to  the  current. 

We  now  investigate  the  effect  of  disorder  on  the 
transmission  probability  through  a  potential  barrier  in 
graphene.  This  can  appear  as  short-range  disorder,  inter- valley 
scattering  and  trigonal  distortion.  In  single  layer  graphene. 
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Figure  3.  Energy  dispersion  and  transmission  probability  for  the  case  of  two  independent  pairs  of  energy  subbands.  Panel  (a)  shows  the 
energy  dispersion  by  taking  into  account  the  coupling  between  the  states  corresponding  to  two  different  photon  occupation  numbers  in  the 
linear  approximation.  In  (b),  we  present  a  schematic  diagram  showing  how  the  transmission  probability  may  be  increased  due  to  the 
presence  of  the  other  energy  subbands.  Panel  (c)  gives  the  transmission  versus  the  incoming  particle  energy  and  the  angle  of  incidence  for 
the  case  of  two  independent  pairs  of  energy  subbands. 
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Figure  4.  Effect  of  disorder  on  the  transmission  probability  using  the  Lorentzian  distribution  model.  In  panel  (a)  we  plot  in  (i)  the 
transmission  probability  as  a  function  of  kx\  for  dressed  electrons  (black  curve)  compared  to  infinite  graphene  (red  curve)  for  head-on 
collision.  The  plot  in  (ii)  shows  a  two-peak  model  distribution  under  the  influence  of  disorder.  In  panel  (b)  we  show  the  transmission 
probability  for  electron  dressed  states  with  A  =  Vq/3  in  the  presence  of  disorder. 


disorder  also  induces  a  metal-insulator  transition  by  creating 
a  dynamical  gap  [22].  Consequently,  this  effect  can  modify 
the  gap  created  by  the  electron-photon  interaction.  In  bilayer 
graphene,  disorder  directly  leads  to  energy  dispersion  with  a 
gap  as  well  as  modifying  the  energy  dispersion  close  to  the 
band  edges.  Additionally,  disorder  may  also  lead  to  localized 
states  inside  a  gap  in  bilayer  graphene  [23,  24]. 

In  this  paper,  we  introduce  disorder  phenomenologically 
through  the  non-conservation  of  the  transverse  electron 
momentum  ky.  As  mentioned  above,  ky  is  conserved  for  both 
Dirac  electrons  in  infinite  graphene  and  states  dressed  by 
photons.  Introducing  the  quantity  F  as  a  measure  for  the 
disorder  [25,  26],  we  model  its  stochastic  distribution  as  a 
Lorentzian,  yielding 


F  f°° 

hl\s(kx?  ky )  -  /  d^y 

J — oo 


tp(kx,  qy) 
(ky-qy)2  +  r2’ 


(34) 


where  to(kx,  qy)  denotes  the  transmission  coefficient  in  the 
absence  of  any  imperfections.  As  long  as  the  disorder  is 
weak  with  F  <$(  1,  different  distributions,  which  give  the 
6-function  in  the  limit  of  F  — »•  0,  will  result  in  almost  equal 
transmission  coefficients.  For  nearly  head-on  collision,  ky  <<C 
kx 2,  the  transmission  coefficient  may  be  obtained  analytically 
using  a  Gaussian  distribution.  The  imperfect  boundary  of  the 
potential  region  can  also  be  the  result  of  some  stochasticity  of 
kX2  to  make  the  effect  stronger.  We  neglect  this  effect  since 
our  goal  is  to  investigate  the  role  played  by  disorder  using  a 
simple  approximation. 


Our  numerical  results  showing  the  effect  due  to  disorder 
are  presented  in  figure  4.  First  we  test  our  numerical  results 
in  (a)  by  applying  the  Lorentzian  transformation  in  (34)  to 
the  complete  transmission  with  T  =  1,  corresponding  to  the 
Klein  paradox  with  head-on  collision  for  infinite  graphene. 
Analytical  integration  clearly  results  in  a  transmission 
amplitude  equal  to  unity.  The  transformed  distribution 
demonstrates  the  precision  of  our  numerical  procedure.  By 
comparing  figure  4(b)  with  figure  1(d),  we  see  a  complete 
suppression  of  the  perfect  transmission  with  T  =  1  along  the 
diagonal  kx\  —  kx 2  by  disorder  along  the  boundary  as  well  as  a 
reduced  range  of  kx 2  for  photon-assisted  perfect  transmission. 

4.  Concluding  remarks 

According  to  recently  published  results  [9],  the  interaction 
between  Dirac  particles  in  graphene  and  circularly  polarized 
light  leads  to  the  formation  of  quantum  electron  dressed  states. 
These  states  are  appreciably  different  from  conventional  Dirac 
electrons  in  ordinary  infinite  graphene.  From  [9],  a  laser 
power  of  100  mW  leads  to  a  gap  of  the  order  of  A  ^  100  meV, 
which  is  required  to  make  the  effect  significant  for  infrared 
light  frequencies  and  room  temperature.  This  enables  possible 
experimental  demonstrations  of  the  described  effects.  We  have 
shown  that  the  electron-photon  interaction  gives  rise  to  states 
of  broken  chiral  symmetry.  The  non-symmetrical  properties  of 
the  states  become  more  significant  when  the  electron-photon 
interaction  is  increased.  In  addition,  there  are  no  dressed  states 
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with  chirality  symmetry.  In  general,  incoming  electrons  or 
holes  passing  unimpeded  through  a  square  potential  barrier 
require  chiral  symmetry.  Under  illumination  from  a  circularly 
polarized  light,  we  can  control  the  degree  of  partially  broken 
chiral  symmetry  in  dressed  states  or  the  degree  of  partially 
perfect  transmission  through  a  potential  barrier  in  a  graphene 
layer. 

In  the  simplest  approximation  when  only  the  two  nearest 
subbands  are  included,  the  model  is  formally  similar  to 
the  so-called  0-3  Hamiltonian  used  to  describe  the  particles 
in  single  layer  graphene  with  parabolic  energy  dispersion, 
giving  non-zero  electron  effective  mass.  We  have  discussed 
the  similarities  as  well  as  the  differences  affecting  the 
wavefunctions  but  not  the  energy  dispersion.  We  have  also 
studied  in  detail  both  numerically  and  analytically  how  the 
transmission  probability  depends  on  the  angle  of  incidence  or 
on  the  incoming  energy  for  the  case  of  head-on  collision. 

By  including  more  of  the  next-nearest  subbands,  the 
tunneling  amplitude  is  modified  in  a  significant  way.  In  the 
approximation  when  two  independent  pairs  of  subbands  are 
included,  we  obtain  an  enhanced  transmission  probability 
when  the  barrier  height  is  close  to  hw o-  This  is  due  to  the 
fact  that  one  of  the  terms  in  the  corresponding  wavefunction 
could  be  perfectly  transmitted.  We  have  analytically  obtained 
the  energy  dispersions  for  this  case.  By  including  more  than 
two  independent  pairs,  this  effect  will  decrease  since  perfect 
transmission  will  occur  only  for  the  two  wavefunction  terms. 
The  effect  is  not  sensitive  to  the  barrier  width,  and,  therefore, 
can  be  considered  as  reminiscent  of  the  Klein  paradox. 

We  have  introduced  disorder  phenomenologically  in  this 
paper  through  non-conservation  of  the  transverse  electron 
momentum  component,  which  is  shown  to  suppress  the 
perfect  transmission  along  the  diagonal  kx\  —  kx 2.  From 
a  physical  point  of  view,  this  disorder  model  could  be 
interpreted  as  arising  from  the  surface  roughness  of  the 
potential  barrier.  The  same  type  of  statistical  distribution 
can  be  applied  to  fluctuations  in  the  barrier  width  which 
will  result  in  modification  of  the  intensity  and  location 
of  the  transmission  peaks.  Consequently,  the  transmission 
maxima  observed  experimentally  will  not  exactly  match  those 
theoretically  predicted  for  clean  samples. 
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